Quantum tunneling between paramagnetic and superconducting states of a 
nanometer-scale superconducting grain placed in a magnetic field. 
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We consider the process of quantum tunneling between the superconducting and paramagnetic 
states of a nanometer-scale superconducting grain placed in a magnetic field. The grain is supposed 
to be coupled via tunneling junction to a normal metallic contact that plays a role of the spin 
reservoir. Using the instanton method we find the probability of the quantum tunneling process and 
express it in terms of the applied magnetic field, order parameter of the superconducting grain and 
conductance of the tunneling junction between the grain and metallic contact. 



Recent advances in manufacturing of small elec- 
tronic devices posed several questions in the theory of 
nanometer-size superconductors related to their prospec- 
tive applications ]jj . One of the key issues is the behavior 
of ultra-small superconductors, of the dimensions much 
less than the coherence length £ in the presence of a mag- 
netic field. Of a special interest from both fundamental 
and practical point of view are the dynamical processes 
of switching between the superconducting and paramag- 
netic states. Small grains can experience spontaneous 
transitions due to quantum fluctuations. In this Letter 
we investigate the process of quantum tunneling between 
the superconducting and paramagnetic states of an ultra- 
small superconductor. 

In geometrically restricted superconductors the Zee- 
man effect suppresses superconductivity at magnetic 
fields exceeding H sp i n = A/V2/UB where A is the order 
parameter |^| Q . The orbital effect in small superconduc- 
tors, which usually plays a major role, was discussed in a 
pioneering work by Larkin Q : In a dirty spherical grain, 
At <C 1, t is the normal excitation relaxation time, the 
orbital effect leads to suppression of superconductivity 
at H orb ~ $o/(r v /D/A), where $ = hc/2e, D is the 
diffusion coefficient, and r is the radius of the grain. 
Tunneling switching effects appear in the paramagnetic 

'l/r <C 1, where I — vpr 



limit H spin /H, 



lorb 



is the mean free path, and d is the mean energy level 
spacing. In the ballistic case (I ~ r) this limit is achieved 
only if A < d corresponding to the strongly fluctuating 
regime || . However in a platelet geometry with magnetic 
field applied along the film H spin / 'H or b ~ yA/d yUb/S, 
where b is the thickness and 5* is the area of the sample; 
thus the paramagnetic limit can be easily achieved along 
with condition A 3> d even in the ballistic case if the 
ratio b/yS is sufficiently small. 

We focus on the case where the Zeeman effect dom- 
inates the orbital one while the order parameter A is 
larger than the mean level spacing d and investigate the 
process of quantum tunneling between superconducting 
and paramagnetic states. Since these two states are char- 
acterized by different values of the total spin, such quan- 



tum tunneling process can occur only in presence of pro- 
cesses allowing for non-conservation of the total spin. We 
consider the grain weakly coupled with a normal metal- 
lic lead (plate) that plays a role of a spin bath. Our 
final result for the probability of the quantum tunneling 
between the superconducting and paramagnetic states is 



P ~ exp [N \n(f38EG/A )' 



(1) 



where the numerical coefficient (3 ss 1.1, the factor N 
is the number of polarized electrons in the paramagnetic 
state of the grain, G is the conductance of the tunneling 
junction between the grain and the metallic lead mea- 
sured in units e 2 /h, A is the order parameter of the 
grain in the superconducting state and NSE is the to- 
tal energy difference between the grain's superconduct- 
ing and paramagnetic states. The factor N is related to 
the applied magnetic filed H and the average density of 
states v as N — 2v^bH with jiB — \e\h/2mc. The en- 
ergy difference per one state 8E is related to the magnetic 
field deviation SH from the magnitude at which the the 
thermodynamic phase transition occurs via SE = fisSH. 
The result (|l|) is valid as long as SE/A <C 1, the general 
case A ~ SE is more complicated and we leave it for 
future study. The meaning of our answer ([!]) can be un- 
derstood as follows. The probability of a single tunneling 
process is (5E/ Aq)G, then the probability of having N of 
them simultaneously in order to transit to paramagnetic 
state is simply [(SE/A )G] N . 

The model. We consider a system of a superconducting 
grain and a metallic plate coupled via weak tunneling 

H = H g + H M + tkk ' tib*** + 4v4a)]- (2) 

k,k' 

where H g and Hm are the Hamiltonians of the grain and 
metallic plate, 4* (■0) and d'(d) are the creation (annihi- 
lation) operators of electrons of the grain and of the metal 
respectively and tkk' is the electron tunneling matrix ele- 
ment between the grain and metal. Electrons of metallic 
plate are described by the frcc-fcrmion Hamiltonian 



Hm — ^ d\,, a C, a k' dk 
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while the grain is described by the BCS model 



fe2T, 



(3) 



k\.k2 



where £/. are the exact eigenvalues (measured with re- 
spect to the chemical potential) of the noninteracting 
Hamiltonian, A is the coupling constant, and the mag- 
netic field h, pointing along the z— axis, is measured in 
the energy units h — \xbH- 

Quantum tunneling. The amplitude of the tunneling 
process between the initial and final states is 



A =< f | T t e 



i J tj H(t) dt 



I >, 



(4) 



where the Hamiltonian H(t) is written in the interac- 
tion representation: the noninteracting part includes the 
metal Hamiltonian Hm and the first (free-fermion) term 
of grain Hamiltonian while the interaction part in- 
cludes the BCS interaction and electron tunneling be- 
tween the grain and metal. We will consider a process of 
quantum tunneling from pure superconducting state (ini- 
tial state) to the paramagnetic state with order parame- 
ter A = 0. The electron tunneling between the grain and 
metal will be treated as a perturbation assuming that the 
tunneling matrix elements tkk/ arc small. The paramag- 
netic state has a nonzero total spin S which is formed by 
the polarized electrons with \£ k \ < £ such that S = 
During the quantum tunneling process the spin of the 
grain increases from zero to S, thus there must be 25* 
electron tunneling processes between the grain and metal 
and first nonzero contribution in expansion of Eq.(|J) in 
tunneling element emerges only in N = 2S order. It 
is clear that the paired states which are destroyed by 
the electron tunneling are those with |£| < £. Expanding 
Eq.(0) it t we have 



A=<f\ 



dt k T t e 



dt ■ 



fe' 



tkk' 



n 

l£fcl<l 

x $L(*fc)4v(*fe) +dl, a (t k )i>ka{tk)] \i > 



(5) 



In the absence of coupling between the grain and metal 
the initial and final states of the system are the products 
of the corresponding initial and final states of the grain 
and metal \i >= \i G > \i M >, 1/ >= |/g > >, 
thus in the leading order in tunneling matrix clement the 
quantum mechanical average of operators d in da) can be 
directly implemented. Since the total spin of the metallic 
plate decreases during the quantum tunneling process the 
only relevant matrix elements are those that correspond 
to creation of electrons with spin down and annihilation 
of electrons with spin up 

<fM\dl l (t)\i M >=e i ^ t , Cfei>0 
</m|<4t(*)Km >=e l ^Tl* ; CfeT <0 . 



The initial state of the metal is the Fermi sea while its 
final state having TV electron-hole excitations is charac- 
terized by the set of excitation energies Q,^ where the 
index k' labels the quantum number of the excitation p. 
Now assuming that tunneling matrix elements are index 
independent tk.k' — t the amplitude can be written as 

A = t"<f G \T t J] fdt k e-*V* Wdt £ e <# W> * 

|ffc|<e Perp(fe) 



^Uh)0{~$ ik)} ) + M(h)0{tir")] \i G > (6) 

where the sum goes over all excitation permutations p(k). 
Since we consider only the case SE < < A we can simplify 
the problem assuming that the energy of a typical excita- 
tion in the metal is much less than A. In this case we can 
neglect the energies (ffi in (^) and using the Heisenberg 
representation we write the amplitude (|^) as 

A = N\t N <f G \ J] jdt k [^(t k )e(-c { k f m ) 

\tk\<£ 

+Mt k )0(cl? {k)} )}\*G>, (7) 

where a specific permutation p(k) was chosen. 

The probability P of the tunneling process is obtained 
by integrating of A* A over all possible quantum states of 
the metal 

N f 

P~ -§ J d( {1} ..d( {N} A*A5(N SE -^2\( {p} \), (8) 



{p(fc)h 



where vm is the density of states of the metal. 



-T t X 

FIG. 1. Dependence of the order parameter A on time 
t. The region r > is the mirror reflection of the "physical" 
region r < 0. The superconducting state corresponds to the 
regions r| > T while the paramagnetic one to |r| < T. 

To find the value of A we shall use the instanton 
method turning to use the Euclidean time t — > —it. Tak- 
ing initial and final times as Tj = — oo, Tt = we present 
the amplitude A as 

A = Nlt N < f G \S{0,-oc)\i G > 
with the evolution operator 



sfo,n)= n 

\tk\<l 
+MTk)e(C { J m )]. 



drkW k] (r k )6(-Cif [m ) 



(9) 
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In principle, the instanton method should be applied to 
the amplitude A directly, but it is more convenient to 
consider the product A* A presenting A* as A* — N\ t N < 
io\S^ (oo, 0)|/g > and writing the product A* A as 

A* A = < i G \ 5t(oo,0) |/ G > < f a \ S(0, -oo)\i G > 
= <i G |5t(o 0) 0)S , (0,-oo)|i > ) 

where it was used that by construction the Euclidean evo- 
lution operator S brings the grain from initial \ic > to 
the final \ fa > state: S \iq >— \ fa > ■ Now the instan- 
ton process has the following structure (see Fig.l.): the 
evolution begins at r = — oo from the superconducting 
state, then at r w — T the system turns into the para- 
magnetic state and stays there till it turns back to the 
superconducting state at r « T. The artificial part of the 
process (r > 0) is the mirror reflection of the "physical" 
process with r < 0. The advantage of this representation 
is that now the new initial (r = — oo) and final (r = oo) 
states are identical and therefore we can use the conve- 
nient functional representation for A* A 

A*A= [N! t N f J Di> el Cdt J] / dr lk dr 2k 

k 

x ^ T (^)^ T (r 2fc )0(-C, { f (fe)} ) 



with the Lagrangian 



(10) 



C = --A*A - V) V|[9 T + a - ha z 

A k 

-A^ T <-A*Vfc|^ T - (11) 

The integration over the fermionic fields in ( |Io| ) should 
be implemented exactly while the assumed integration 
over A(t) will be done with the saddle point accuracy. 
Now taking integrals in d§) for the probability we obtain 



p _ e Nln(N SEt 2 v M ) 



k 

xel CdT [^ k \{Ti k )^{T 2k ) +ipli(nk)ipki{T2k)], 
and finally integrating over the fermionic fields we obtain 

In P = N ln(N 5E t 2 v M ) + Tr ln t 9r + H ^ 

k 

-- f drA*(r)A(r)+^lnZ fc , 
9 J j. 



(12) 



where Z k = Tr J dndr 2 G lk {ri,T 2 ), 

A(r) 



Wfc(r) = 



A*(r) 



(13) 



and the matrix Green function 



Gik(n,T 2 ) 



Gi k (Ti,T 2 ) Fifc(ri,r 2 ) 
FUn,r 2 ) G lk {r u r 2 ) 



is defined by the equation 

[dri + Hfc(ri)] G lk {r U T 2 ) = 5(n - r 2 ). 



(14) 



(15) 



Instanton equations. To find the instanton equations 
we take the functional derivative of Eq. ( |l2| ) with respect 
to A* obtaining 



A(r) 



k 



hk(r) + f 2k (r), 



(16) 



where the Green function fx (r) = F± (r, r) emerges form 
the variation derivative of Trln[9 r + H k ] in ([l2]) and f 2 
emerges from the functional derivative of the Green func- 
tion G\ 



/ 2 (r) 



SA*(t 



■ J dridr 2 Tr G lk [n , r 2 ) 



(17) 



Combining Eq.(|T^) with the same equation written in 
the transposed form one can find the equation that de- 
fines the function f x in terms of only diagonal Green 
functions 



d T 9ik(r) + [Hok(r),gik{T)] = 0, 



where 



9lk(r) 



9lk(r) /ijfe(r) 
f!k(r) 9ik(r) 



Gifc(r,r), 



(18) 



(19) 



and H is given by Eq.(|13|) with h = 0. Writing Eq.fll* 
in components we get 



drgik + A/ x + fc - A*/l* = 0, 
SrAk + - 2A.9u. = 0, 

-9 T 4 + 2U\ k - 2A*.g lfc = 0, 
d T s z \k = 0, 



(20) 
(21) 



where the variables gi = \g X k - ffifc]/2, s 2 ifc = -[gik + 
gxk]/2 were introduced instead of components g\ and 
g~\, Eqs .(|2C| ) are very similar to the well known Eilen- 



berger |7j equations and posses the same invariant g\ k 
fl k flk = const. Now we turn to the function f 2 (r) : Us- 
ing the definition of the Green function ( |l5|) one can take 
the variational derivative in Eq.(^) obtaining 

hk{r) = -^l/fMsiW+sfM/^T)] (22) 

where gl, fl are the components of the matrix Green 
function 



9&t) 



9i(r) flC 



fkir) gi(r) 



dT 2 G lk (T,T 2 ) (23) 



and f I 1 , gj. 1 a.re the components of the matrix Green func- 
tion g^{j) — J driGik(Ti, t). Equations that determine 
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the functions g 1 and g 11 can be easily found by integrat- 
ing Eq.(|l5|) over t 2 and transposed Eq.(|i~5|) over t\ 



d T g I k (T)+n k (T)g I k (r) = l 



d r g I k I (T)+g 1 k J (r)n k (r) = l. 



ji 



(24) 



Initial and final conditions. Initial and final condi- 
tions should be formulated for the Green function g k (r) 
which is defined with respect to the evolution operator 
fits, 



~Ot,8 / \ 



< 



i G |5t(oo, 0) 5(0, T)faffis(T, -co) \ia > 



<i G \&(oo,0)S(0,-oo)\i a > 



where a negative r was chosen for concreteness. To find 
the function g k (r) in the functional representation we 
add the source term f c?T/i" /3 (r)?/' Q .(r)'0^(r) where is 
the Nambu spinor to the Lagrangian (|ll]) and take the 
variation derivative of (|l2] ) with respect to /Lt^(r) ob- 
taining 

fffc(-r) = Sifc(r) +g 2 fc(-r), 
where gu- is defined by Eq.([19|) and the function 



hk{r) = 



g2k(r) f 2 k(r) 
fiki T ) 92k(r) 



is related with g^and g 11 by 

-glgV/Zk. 



92k 



(25) 



In the superconducting phase (|r| ^> |T|) the function 
<7fc should coincide with the equilibrium superconducting 
Green functions 

g k = ^ k /2E k , f k = A/2E k , ft = ,s z = (26) 



where £ fe = v£| + A 2 , g fc = [ff* - <7fe]A s zk = -[g k + 
g~k]/2. The function g is directly related with electron 
density on the level k by n k = l — 2g k and the function s z k 
is the z-component of the spin on the level fc. Analogously, 
in the paramagnetic phase the Green function g is 



9k = fk = fk = 
g k = sign^ fc /2, 



0. 

Szk 



Szk 
= fk 



1/2 

"ti- 



ter |ffc| 
for \£k\ 



(27) 
(28) 



In the absence of tunneling the physical Green function 
g would coincide with gi that obeys equations ( p0|j2l| ) 
conserving the quasiparticle spin. Therefore the function 
gi in the paramagnetic state obeys the boundary condi- 
tion fl2S| ) for any ^ while in the superconducting state it 
obeys the boundary conditions (pGj), so that §2 — > in 
the superconducting phase. 

Numerical solution. Solution of Eqs.(pO|) satisfying the 
necessary boundary conditions for a given configuration 
of the order parameter A(r) can be easily found numeri- 
cally. To find the components of the function 52 one first 




FIG. 2. The normalized order parameter A(t)/Ao (dashed 
line) and the normalized total spin 2S(r)/N of the grain as 
functions of on r' = r — T at the boundary r ~ T between 
the superconducting (left) and paramagnetic (right) states. 



needs to solve Eqs.(f24|) numerically and then find the 
components §2 according to Eq. (|25|) . Knowing the func- 
tions <7i and 32 for a given configuration A(t) one can 
find the self-consistent configuration of A(r) satisfying 
Eq. (|l6|) . This self-consistent solution A(r) along with the 
total spin of the grain S(t) at the right instanton bound- 
ary is presented on Fig. 2 (solution at the left boundary 
can be obtained as the mirror reflection of that on the 
right one). Substituting the functions A(r) and g 1 \g TI 
into Eq.([T^) and calculating the term Trln[9 T + Hk] by 
the method described in |g| one obtains the result (|l|). 

Conclusions. In conclusion, we have found the proba- 
bility of the quantum transition between superconduct- 
ing and paramagnetic states of a nanometer size super- 
conducting grain weakly coupled to a normal metallic 
plate. Our result (|l|) obtained formally in the limit of 
T — > is general and remain valid at finite temper- 
atures as long as T <C SE. At higher temperatures, 
SE <C T <C A, the characteristic energy 5E in the ex- 
ponent of Eq.(|l|) has to be substituted by temperature 
T. However, finding numerical coefficient (3 in this case 
requires a more advanced study. 

We would like to thank Y. Bazaliy, Ya. M. Blanter, 
E.M. Chudnovsky, D. Feldman, Y.M. Galperin, A.E. 
Koshelev, A.I. Larkin, K. Matveev, A. Mel'nikov, V.L. 
Pokrovsky, and R. Ramazashvily for useful discussions. 
This work was supported by the U.S. Department of En- 
ergy, Office of Science under contract No. W-31-109- 
ENG-38. 



[1] C.T. Black, D.C. Ralph, and M. Tinkham 76 (1996) 688 

[2] A. M. Clogston, Phys. Rev. Lett., 5 464 (1960) 

[3] R. Meservey, P.M. Tedrow, and P. Fulde, Phys. Rev. Lett. 

25 (1970) 1270 
[4] R. Meservey, P.M. Tedrow, Phys. Rep. 238 (1994) 173 
[5] A.I. Larkin, Zhurnal Eksperimental'noi i Teoreticheskoi 

Fiziki 48 232 (1965); (Sov.Phys. JETP 21, 153 (1965)). 
[6] P. W. Anderson, J. Phys. Chem. Solids 11 (1959) 28 
[7] G. Eilenberger, Z. Phys. 214, 195 (1968) 
[8] A.V. Lopatin and L.B. Ioffe, Phys. Rev. B 6412 (1999) 



4 



